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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

feet 

acres 

cubic  yards 
tons  (2000  lb,  mass) 
feet  per  second 
square  feet  per  second 

feet  per  second  per  second 




0.301+8 
1+01+6.856 
0.761+551+9 
907. 18U7 
0.301+8 
0.09290301+ 

0.30U8 


To  Obtain 

metres 

square  metres 

cubic  metres 

kilograms 

metres  per  second 

square  metres  per 
second 

metres  per  second 
per  second 
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LONG  BEACH  HARBOR  NUMERICAL  ANALYSIS  OF  HARBOR  OSCILLATIONS 


ALTERNATE  PLANS  FOR  PIER  J COMPLETION  AND  TANKER 
TERMINAL  PROJECT  (62-  AND  82-FT  DEPTHS) 


PART  I : INTRODUCTION 


1.  The  objective  of  this  study  was  to  investigate  by  use  of  a 
numerical  model  the  response  to  long-period  wave  excitation  of  six  pro- 
posed alternatives  for  modifications  to  Pier  J in  Long  Beach  Harbor  for 
oil  tanker  berthing  and  general  cargo  facilities.  Wave-height  amplifi- 
cation factors  and  normalized  maximum  current  velocities  were  plotted 
versus  wave  period  at  the  proposed  oil  tanker  terminals  to  ascertain 
whether  or  not  significant  harbor  oscillations  might  occur  which  would 
pose  problems  to  ship  berthing. 


2.  The  tanker  terminal  configurations  (plans)  were  considered  in 
the  study;  plan  1 is  shown  in  Plate  1.  Both  300-  and  600-ft*  openings 
at  the  east  end  of  the  basin  and  62-  and  82- ft  depths  were  tested.  The 
basin  geometry  was  slightly  different  for  the  82-ft  depth  since  there 
was  a cutback  with  dimensions  of  approximately  l60  by  1550  ft  at  the 
east  end  of  the  basin  (Plate  2).  Plan  2 (Plate  3)  was  similar  to  plan  1 
except  that  the  proposed  landfill  addition  to  Pier  J was  reduced  to 
approximately  50  acres  with  an  open-water  area  to  the  east  for  mooring 
two  250,000-dwt  tankers.  Water  depths  of  62  and  82  ft  were  tested. 

Waves  incident  from  a direction  normal  to  the  main  breakwater  face  with 
periods  from  30  sec  to  300-3^0  sec  were  considered  for  both 


i H ftt 


configurations.  Wave  amplitude  and  current  velocities  were  calculated 
every  5 sec  for  the  period  range.  Resonant  peaks  also  were  defined  by 
considering  incident  wave  periods  in  increments  as  low  as  0.25  sec. 


PART  II:  NUMERICAL  MODEL 


3.  The  response  of  the  tanker  terminal  basin  to  long  wave  excita- 
tion was  determined  by  using  a hybrid  finite  element  numerical  model  de- 
veloped recently  by  Chen  and  Mei*  at  the  Massachusetts  Institute  of  Tech- 
nology. The  model  solves  the  following  generalized  Helmholtz  equation: 

2 

V • [h(x,y)V<|>(x,y)]  + <p(x,y)  = 0 (l) 

o 

where  <j>(x,y)**  is  the  velocity  potential  defined  by  u(x,y)  = -V(f>(x,y), 
with  u(x,y)  being  a two-dimensional  velocity  vector  and  w an  angular 
frequency.  Equation  1 governs  small  amplitude  undamped  oscillations  of 
water  in  a basin  of  arbitrary  shape  and  variable  depth  forced  by  periodic 
long  waves.  It  has  been  further  assumed  that  the  flow  is  irrotational. 

U.  The  boundary  condition  along  the  shoreline  and  along  the  de- 
tached breakwater  surrounding  the  three  proposed  berths  is  that  the 
normal  component  of  the  velocity  be  equal  to  zero.  Therefore,  the 
breakwater  is  considered  as  a solid  barrier.  No  special  boundary  con- 
ditions are  made  to  account  for  the  storage  tankers  in  plan  2 or  any 
vessels  utilizing  the  tanker  terminal  basin,  since  for  long-period  waves 
of  small  amplitude  ship  heave  equals  incident  wave  height.  Therefore, 
the  vessels  do  not  significantly  alter  the  oscillation  characteristics 
of  the  basin. 

5.  The  Helmholtz  equation: 

2 w2 

V <|>(x,y)  + — <f>(x,y ) = 0 (2) 

is  the  governing  equation  for  a constant-depth  ocean  region  outside 
the  basin. 

* H.  S.  Chen  and  C.  C.  Mei,  "Oscillations  and  Wave  Forces  in  an  Off- 
shore Harbor  (Applications  of  the  Hybrid  Finite  Element  Method  to 
Water-Wave  Scattering),"  Report  No.  190,  197^,  Massachusetts  Institute 
of  Technology,  Cambridge,  Mass. 

**  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 


6 


6.  For  a harbor  in  a semi -infinite  ocean  with  a straight  coast- 
line there  is  an  incident,  reflected,  and  scattered  wave.  The  scattered 


wave  has  a velocity  potential 


where  a are  unknown  coefficients  and  H (kr)  are  Hankel  functions 
n n 

of  the  first  kind  of  order  n . 

7-  <b  satisfies  the  radiation  condition  that  the  scattered  wave 
must  behave  as  an  outgoing  wave  at  infinity.  This  condition  is  known 
as  the  Sommerfeld  radiation  condition  and  may  be  expressed  mathemat- 
ically as  follows: 


8.  Chen  and  Mei  used  a calculus  of  variations  approach  and  ob' 
tained  a Euler-Lagrange  formulation  of  the  boundary  value  problem.  1 
following  functional  with  the  property  that  it  is  stationary  with 
respect  to  arbitrary  first  variations  of  <j>(x,y)  was  constructed  by 
Chen  and  Mei : 


F(<j>)  = ff  1/2  (h(V<f>) 


A = region  inside  the  harbor 
§ = line  integral 
1 = far  field  velocity  potential 


$ = velocity  potential  of  the  incident  wave 

n = unit  normal  vector  outward  from  region  A 

SL 

a = boundary  of  region  A 

<(>  = total  velocity  potential  evaluated  on  boundary  a 

8. 

9.  Proof  was  given  by  Chen  and  Mei  that  the  stationarity  of  this 
functional  is  equivalent  to  the  original  boundary  value  problem. 

10.  The  integral  equation  obtained  from  extremizing  the  func- 
tional is  solved  by  utilizing  the  finite  element  method.  This  method  is 
a technique  of  numerical  approximation  that  involves  dividing  a domain 
into  a number  of  nonoverlapping  subdomains  which  are  called  elements. 

11.  The  solution  of  the  problem  is  approximated  within  each  ele- 
ment by  suitable  interpolation  functions  in  terms  of  a finite  number 
of  unknown  parameters.  These  unknown  parameters  are  the  values  of  the 
field  variable  <J>(x,y)  at  a finite  number  of  points  which  are  called 
nodes.  The  relations  for  individual  elements  are  combined  into  a 
system  of  equations  for  all  unknown  parameters. 

12.  In  the  region  outside  the  basin,  the  velocity  potentials  are 
solved  analytically  in  terms  of  unknown  coefficients.  The  region  is 
considered  a single  element  with  an  "interpolation  function"  given  by 
Equation  3.  The  infinite  series  is  terminated  at  some  finite  value  such 
that  the  addition  of  further  terms  does  not  significantly  influence  the 
calculated  values  of  <j>(x,y)  . The  resulting  equation  is  combined  with 
the  system  of  equations  for  unknown  parameters  at  nodal  points  within 
the  basin  and  this  complete  system  is  solved  using  Gaussian  elimination 
matrix  methods. 


13.  n(x,y)  is  related  to  4>(x,y)  through  the  linearized  dynamic 
free  surface  boundary  condition 


15.  The  hybrid  finite  element  method  (so  named  by  Chen  and  Mei 
because  the  method  involves  the  combination  of  analytical  and  finite  ele^ 
ment  numerical  solutions)  is  a steady-state  solution  of  the  boundary 
value  problem.  The  response  of  a harbor  to  an  arbitrary  forcing  func- 
tion can  be  easily  determined  within  the  framework  of  a linearized 
theory. 


PART  III : RESULTS 


16.  Plate  4 shows  the  locations  of  the  tanker  terminals  and 

Plates  5-7  show  the  finite  element  grids  for  STFP  2 (600)  with  62- ft 
depth,  STFP  2 (600)  with  82- ft  depth,  and  STFP  3.  The  grids  for  STFP 
(300)  with  62-  and  82- ft  depths  are  not  shown  hut  are  similar  to  the  grids 
for  STFP  (600)  except  for  the  breakwater  opening.  Plates  8-13  are  plots 
of  normalized  maximum  current  velocity  versus  wave  period  at  each  tanker 
terminal  for  the  six  modifications.  The  current  velocity  multiplied  by 
the  incident  wave  amplitude  in  feet  gives  velocity  in  units  of  feet  per 
second.  The  velocities  have  no  vertical  component  or  variation  since  only 
long  waves  are  considered  in  this  study.  Also,  the  maximum  velocity  over 
one  wave  period  is  plotted.  Plates  14-19  are  plots  of  the  wave-height 
amplification  factor  versus  wave  period  in  seconds  at  each  tanker  termi- 
nal. The  wave-height  amplification  factor  is  defined  at  a point  inside 
the  basin  as  the  wave  height  at  the  point  divided  by  twice  the  incident 
wave  height.  This  definition  of  amplification  factor  is  traditional  and 

is  a result  of  the  fact  that  the  standing  wave  height  for  a straight 

coast  with  no  harbor  would  be  twice  the  incident  wave  height  due  to  the 
superposition  of  the  incident  and  reflected  waves. 

17.  Plan  STFP  2 (300)  and  plan  STFP  2 (600)  for  the  62- ft  depth 

have  normalized  maximum  current  velocity  (UMCV)  resonant  peaks  of  57  and 
280  sec  and  57*5  and  260  sec,  respectively.  The  NMCV  patterns  in 

Plates  20  and  22  for  the  peaks  of  approximately  1 min  are  very  similar 

for  the  two  plans  with  the  velocities  somewhat  larger  for  plan  STFP  2 

(600).  The  NMCV  patterns  in  Plates  21  and  23  for  the  peaks  of  approxi- 

mately 4-1/2  min  are  virtually  identical. 

18.  Plan  STFP  3 has  NMCV  resonant  peaks  of  55,  55*5,  106,  and 

350  sec.  The  NMCV  peaks  of  approximately  1 min  are  substantially  larger 
than  similar  peaks  for  plan  STFP  2 (300)  and  STFP  2 (600).  The  NMCV  pat- 
terns in  Plates  24  and  25  are  similar  to  the  patterns  in  Plates  20  and  22 

with  high  velocities  occurring  near  corners  of  the  breakwater.  Large  ve- 
locities also  occur  in  the  area  of  the  storage  tankers.  The  106-sec 
NMCV  peak  is  substantially  larger  for  plan  STFP  3 than  for  the  other  two 


plans;  Plate  26  shows  the  NMCV  pattern  for  this  peak.  Velocities  are 
very  large  at  terminals  2 and  3 as  a result  of  the  antinode  forming  in 
the  inside  corner  of  the  breakwater  (see  Plate  Uo  for  the  contours  of 
wave-height  amplification).  The  350-sec  NMCV  peak  at  terminal  1 is  simi- 
lar in  magnitude  to  the  U-l/2-min  peaks  at  terminal  1 for  the  other  two 
plans.  However,  the  NMCV  peaks  at  terminals  2 and  3 are  less  than 
corresponding  peaks  of  plans  STFP  (300)  and  STFP  (600).  Plate  27  shows 
the  NMCV  pattern  for  the  350-sec  peak. 

19*  Wave-height  amplification  factors  are  quite  similar  for 
plans  STFP  (300)  and  STFP  (600)  (Plates  l4-l6).  Contours  of  wave-height 
amplification  presented  in  Plates  3U-37  are  also  very  similar.  Antinodes 
for  the  1-min  peak  are  located  at  the  two  inside  corners  of  the  detached 
breakwater  and  along  the  shore  near  terminal  3.  A single  antinode  in 
the  center  of  the  basin  forms  for  the  4-1/2-min  peak. 

20.  The  1-min  peaks  for  STFP  3 are  substantially  larger  than 
corresponding  peaks  for  the  other  two  plans.  Antinodes  form  in  the  in- 
side corners  of  the  breakwater,  near  terminal  3,  and  along  the  corner 
of  the  landfill  (Plates  38  and  39)-  Contours  of  wave-height  amplifica- 
tion for  the  350-sec  peak  of  STFP  3 (Plate  4l)  are  similar  to  the  con- 
tours for  plans  STFP  2 (300)  and  STFP  2 (600). 

21.  Plans  STFP  2 (300)  and  STFP  2 (600)  for  the  82- ft  depth  have 
NMCV  resonant  peaks  (Plates  28-31)  and  wave-height  amplification  peaks 
(Plates  h2-h'jj)  similar  to  the  peaks  for  the  62-ft  depth  except  that  the 
peaks  for  the  82- ft  depth  are  somewhat  larger.  This  difference  appears 
to  be  caused  by  the  geometry  change  of  the  cutback  with  dimensions  of 
approximately  l60  by  1550  ft  at  the  east  end  of  the  basin.  For  ex- 
ample, a comparison  of  Plates  3U  and  k2,  which  show  contours  of  wave- 
height  amplification  for  the  1-min  peaks  for  STFP  2 (300),  reveals  that 
the  cutback  has  increased  the  height  of  the  antinode  along  the  landfill 
near  terminal  3.  The  antinode  in  the  corner  of  the  detached  breakwater 
has  also  increased  for  the  cutback  conditions. 

22.  The  contours  of  wave-height  amplification  (Plates  U6  and  U7) 
and  NMCV  patterns  (Plates  32  and  33)  are  almost  identical  for  the  62- 
and  82- ft  depths  of  STFP  3;  the  amplification  factors  and  velocities 


are  larger  for  the  62-ft  depth,  however.  This  difference  may  result 
from  greater  reflection  of  energy  for  the  82-ft  depth  as  waves  propagate 
from  the  60-ft  depths  outside  the  basin  to  the  82-ft  depths  inside.  The 
same  effect  might  contribute  to  the  lower  amplification  factors  and 
velocities  for  the  82-ft  depth  of  plans  STFP  2 (300)  and  STFP  2 (600). 

23.  In  summary,  larger  1-  and  2-min  resonant  oscillations  occur 
for  plan  STFP  3 than  for  plans  STFP  2 (300)  and  STFP  2 (600).  The  U- 
to  6-min  oscillations  are  somewhat  larger  for  plans  STFP  2 (300)  and 
STFP  2 (600)  than  for  plan  STFP  3.  The  oscillation  characteristics  of 
plans  STFP  2 (300)  and  STFP  2 (600)  are  very  similar.  Oscillations  are 
also  somewhat  greater  for  all  plans  for  the  62-ft  depth  than  for  the 
82-ft  depth  with  the  l60-ft  cutback. 


PART  IV:  CONCLUSIONS 


2h.  The  resonant  peaks  of  the  normalized  maximum  current  veloc- 
ity and  the  wave-height  amplification  factor  do  not  appear  to  preclude 
safe  berthing  at  the  terminal  areas  given  proper  mooring  practices  and 
adequate  mooring  equipment.  The  resonant  response  peaks  for  incident 
waves  with  periods  of  approximately  1 to  2 min  are  significantly  larger 
for  the  62- ft  depth  STFP  3 plan  than  for  the  62-ft  depth  STFP  2 plans. 
Whether  or  not  incident  waves  in  this  period  range  excite  significant 
ship  motion  depends  upon  ship  characteristics  (length,  width,  and  draft) 
and  mooring  conditions  (type  of  lines,  amount  of  slack,  etc.). 

25.  The  incident  wave  spectrum  and  ship  response  versus  wave 
period  curves  must  be  known  in  addition  to  the  harbor  response  curves  of 
Plates  8-19  to  definitively  determine  whether  or  not  ship  surge  problems 
resulting  from  harbor  resonance  will  occur  in  the  tanker  terminal  basin. 

26.  Until  detailed  information  on  moored  ship  response  as  a func- 
tion of  incident  wave  amplitude  and  frequency  is  obtained,  definitive 
conclusions  on  the  precise  amount  of  ship  motion  cannot  be  made.  How- 
ever, the  oscillations  of  the  tanker  terminal  basin  are  significantly 
smaller  than  typical  oscillations  noted  throughout  the  Los  Angeles  and 
Long  Beach  Harbors  complex  in  the  numerical  study  performed  for  the 
entire  harbor  complex.* 


* J.  R.  Houston,  "Long  Beach  Harbor  Numerical  Analysis  of  Harbor 

Oscillations;  Report  1,  Existing  Conditions  and  Proposed  Improvements" 
(in  preparation),  U.  S.  Army  Engineer  Waterways  Experiment  Station, 

CE,  Vicksburg,  Miss. 
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In  accordance  with  19  70-2-3,  paragraph  6c(l)(h), 
dated  15  February  1973,  a facsimile  catalog  card 
In  Library  of  Congress  format  la  reproduced  below. 
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Long  Beach  Harbor  numerical  analysis  of  harbor 
oscillations;  Report  3:  Alternate  plans  for  Pier  J 
completion  and  tanker  terminal  project  (62-  and  82- 
ft  depths),  by  James  R.  Houston.  Vicksburg,  U.  S. 

Army  Engineer  Waterways  Experiment  Station,  1976. 
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Waterways  Experiment  Station.  Miscellaneous  paper 
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